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Unconventional dc transport in Rashba electron gases 
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We discuss the transport properties of a disordered two-dimensional electron gas with strong 
Rashba spin-orbit coupling. We show that in the high-density regime where the Fermi energy 
overcomes the energy associated with spin-orbit coupling, dc transport is accurately described by a 
standard Drude’s law, due to a non-trivial compensation between the suppression of back-scattering 
and the relativistic correction to the quasi-particle velocity. On the contrary, when the system enters 
the opposite dominant spin-orbit regime, Drude’s paradigm breaks down and the dc conductivity 
becomes strongly sensitive to the spin-orbit coupling strength, providing a suitable tool to test the 
entanglement between spin and charge degrees of freedom in these systems. 
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Spin-orbit (SO) coupling is a fundamental ingredient 
in spintronics[l] as it provides an advantageous locking 
between spin and electron orbital momentum. Recently, 
intense research efforts [2] have been devoted to two- 
dimensional (2D) materials with broken inversion sym¬ 
metry, where the SO strength, parametrized by a char¬ 
acteristic energy scale, Eq, can be tuned by means of ex¬ 
ternal conditions (electric fields, gating, doping, pressure, 
strain, ...). In most of these systems (ex.: surface alloys 
[3-9], layered bismuth tellurohalides [10-16], HgTe quan¬ 
tum wells [17], interfaces between complex oxides [18- 
32]) the total charge carrier density n can be tuned down 
to very small concentrations, implying very small Fermi 
energies Ep. Although the higlr-density (HD) regime 
Ep > E 0 has been widely investigated[2, 33-38], rela¬ 
tively less attention has been paid to the opposite regime 
of dominant SO (DSO), E 0 > Ep. 

In this Letter we provide a detailed investigation of 
the dc conductivity of a 2D electron gas (2DEG) with 
Rashba[39] SO coupling in the different density regimes. 
Using a Boltzmann approach and a fully quantum anal¬ 
ysis based on Kubo formula, we show that in the high- 
density regime Ep > Eq, dc transport is independent of 
the SO strength, and the dc conductivity <Jdc of electrons 
having effective mass m and scattering time r 0 follows 
the conventional Drude law for 2DEGs, 
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that results from a non-trivial cancellation of the SO cou¬ 


pling effects on the quasiparticle velocity and transport 
scattering time. Remarkably, as soon as the system en¬ 
ters the DSO regime E 0 > Ep, Drude’s paradigm (1) 
breaks down and the dc conductivity accurately follows 
the analytical formula: 
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where no = 2mEo/(TrH 2 ) is the density at Ep = Eq. 
In contrast to the linear dependence of ad c on the charge 


density found in the HD regime, n > no, Eq. (2) predicts 
an unconventional non-linear behavior of er^c with n, that 
is controlled by the SO interaction encoded in Hq. The 
relevance of this result is twofold: demonstrating that dc 
transport is strongly sensitive to Rashba SO coupling, 
not only it suggests that SO coupling could be measured 
in a transport experiment but also, what is more impor¬ 
tant for applications, it points at the possibility of tuning 
the conductivity of a 2DEG by tuning the SO coupling 
strength via external gates. 

A disordered Rashba 2DEG confined to the (x,y)~ 
plane is described by the following Hamiltonian 

H =Jidr 'U(r) + a z • (p x <x) + U imp (r) ’P(r), 

( 3 ) 

where Ui mp (r) is the disorder potential, a is the SO cou¬ 
pling, a is the vector of Pauli matrices, 'L(r) and (r) 
are spinor fields which respectively create and destroy 
electrons at position r and a = a/|a|. Here we limit 
ourselves to the simplest case of Gaussian random dis¬ 
order with “white noise” correlations, namely we set 
(Hi mp (r)y imp (r , ))imp = nivf mp S (r —r 7 ) where n imp and rq 
denote respectively the scattering strength and the im¬ 
purity density. In the absence of impurities, H gives an 
electronic spectrum, depicted in Fig.l, consisting of two 
bands with dispersion E* — E 0 = (p + spo) 2 /(2m) — E 0 , 
where po = ma, Eq = ma 2 /2 and s = ±1 denotes the 
eigenvalue of the helicity operator, S, defined as usual as 
S = z • (p x a). In the following we measure the Fermi 
energy, Ep, from the lower band edge, so that Ep = E 0 
corresponds to the “Dirac point”, p = 0, and the HD and 
DSO regimes are realised respectively for Ep > Eq and 
Ep < Eq. 

As already noted in [40], the onset of the DSO regime 
leads to a rather sharp change in density of states (DOS) 
at the Fermi level. Indeed, while for Ep > Eq the DOS 
is constant, N(Ep) = N 0 = m/(ph 2 ), for Ep < Eq it 
displays a van Hove singularity, N(Ep) = Nq^/Eq/E p 
leading to a modification of the dependence of the density 
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Figure 1: (Color online) (a): dispersion of the Rashba model. 
The solid/dashed contours and the red/blue colours denote 
respectively s and 77 . (b): Examples of allowed and forbid¬ 
den back-scattering processes. Thick (red and blue) and thin 
green arrows denote the quasi-particle velocity v ps and the 
spin (<r)ps of each state, respectively. 


on the Fermi energy [40]: 

f Nq(Ef + Eq) Ep > Eq 
[ 2N 0 \/EfEo Ep < Eq. 


(4) 


In the presence of static diluted disorder, the singular 
behavior of the DOS reflects directly on the quasiparticles 
lifetime that, using Fermi golden rule (see e.g. Ref.[4]), 
can be cast as t{E p ) = V E p / a / Qpf s ,] _1 , where V is the 
2D volume of the sample and 

Q P *„> = m i v imp0- + SS 'P -p') S ( E p- E p '). ( 5 ) 

is the scattering kernel. Using Eqs. (4-5) one can easily 
show that the quasi-particle scattering time scales lin¬ 
early with the density in the DSO regime, i.e. 


r = t{E f ) = 


tq Ep > E 0 

tq^/Ep/Eq = r 0 (n/no) Ep < E 0l 


( 6 ) 

where r 0 = H 2 / (mniV 2 mp ) denotes the quasi-particle scat¬ 
tering time in the absence of SO. 

To explain the behavior of the conductivity across the 
different regimes, we start by recalling the definition of 
the velocity operator in the helicity basis, 


[v] ss / = v p As' - ias(l - 5 ss ')i. 


pi 


(7) 


where v ps = V p U* = p(p/m -(- sa) denotes the quasi¬ 
particle velocity and t p is defined as t p = {p y /p, —p x /p}- 
An important thing to underline here, general for any 
chiral system, is that in the presence of SO coupling also 
the velocity acquires a spin structure, which has a deep 
impact on the transport properties. As we discuss in 
more details below, to a first approximation the conduc¬ 
tivity can be described within a standard semiclassical 


Boltzmann approach that only keeps the quasiparticle 
current, arising from the diagonal components of the ve¬ 
locity operator (7). In the relaxation time approximation 
at T = 0, (Jdc can be then estimated as 


a dc ^ o-dc = |y J2 ps 5 ( Ef ~ E p)Wp>K ( 8 ) 

where the transport scattering times T p r s satisfy the fol¬ 
lowing equations [43] 




(9) 


Using explicitly the definition of (5), we introduce 
the transport helicity index p = s(v ps ■ p) = z ■ (v ps x 
(<j) ps ) = ± 1 , that accounts for the reciprocal orientation 
of spin and velocity, and we recast Eq.(9) as follows: 
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where Vf = \/2mEp and p v = \mvF~ppo\ are the Fermi 
momenta on the inner and outer Fermi surfaces. The 
above equation suggests that the index 77 can be used to 
efficiently classify the states at the Fermi level across the 
different regimes. In particular, as illustrated in Fig.la, 
where the value of 77 is indicated by the red/blue col¬ 
ors of the surface, at Ep > Eq p simply coincides with 
s, on the contrary, at Ep < E 0 p allows to distinguish 
between the two Fermi circles that have the same value 
of s but antiparallel quasi-particle velocities. Using this 
classification af c can be cast as 
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where the transport scattering times r* r are given by the 
solution of Eqs. (10) 


= T Prj/PF, (12) 

with pp = 1/2^ i.e. Pf = itivf for Ep > Eq and 
Pf = Po f° r e f < E o- As one can easily check, Eqs.(11) 
and (12) yield Drude result (1) at Ep > E 0 , and Eq.(2) 
in the DSO regime. 

The physical relevance of 77 , as compared to the stan¬ 
dard helicity, s, is evident in Fig.lb, where we show 
that spin conservation forbids back-scattering between 
states having the same value of 77 . What is more impor¬ 
tant, in the DSO regime the only allowed back-scattering 
processes reverse the sign of the quasi-particle velocity 
without changing the direction of momentum. These ef¬ 
fects determine the density dependence of the scattering 
times (a) and of the conductivities (b) of the the major¬ 
ity (v = — 1 ) and minority (77 = + 1 ) carriers shown in 
Fig.2. As one can see in Fig.2b, transport is in general 
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Figure 2: (Color online) (a) Density dependence of the ratio 
t^/t for the majority (rl r , solid blue line) and minority (t+, 
solid red line) carriers, and of the ratio t/to (dashed line), 
(b) Density dependence of Boltzmann dc conductivity a R c , 
and contribution of the two types of carriers, in units of 
oVio = n 0 e 2 To/m. 


dominated by the majority carriers that, due to the sup¬ 
pression of backscattering have also the larger transport 
scattering time, > r > Let us focus on the trans¬ 
port properties of majority and minority carriers across 
the different regimes. In the HD regime r is a constant 
and rji r /r increases as the density decreases due to the 
shrinking of the inner Fermi circle. At n = n 0 , where 
only states with r] = — 1 are present, back-scattering is 
completely suppressed and one recovers t\ L r = 2 r, like e.g. 
in graphene [42]. However, differently from graphene, as 
long as n > no a compensation between SO effects on the 
velocity and on the transport scattering times of the two 
types of carriers restores the usual Drude conductivity, 
even for n very close to no- This result is non trivial: 
indeed, setting naively r* r = r in Eq.( 8 ), would lead to 

[43] (7 — (n — no/2)ro/m i.e. a < dDrude even at n > no 

[44] . On the other hand, as the system enters the DSO 
regime r starts to decrease linearly as predicted by Eq. 
( 6 ) and backscattering processes for the majority carriers 
are progressively restored. Both these effects quench r*L 1 
as ri < no, see Fig. 2a, leading to an overall sublinear 
behavior of the conductivity, see Fig. 2b. 

A deeper insight on dc transport comes from the gen¬ 
eralisation of Boltzmann approach to fully include quan¬ 
tum effects. To this end, we use Kubo linear response 
theory [ 1 ]. We start by noticing that, within the self- 
consistent Born approximation (SCBA), the retarded 
Green’s function is diagonal in the helicity basis and it is 
given by the following matrix [1, 48], 


[G fl (p,w)] ss , = g?<j>,w)8 sa >, (13) 


where g s (p , u) = [w — E p + Ep — E fl (w)] 1 denotes the 
Green function of electrons with helicity s. The self¬ 
energy E r (uj) = niV? mp /(2V) Ep, s fff (P> w ) is spin and 
momentum independent [1], At zero frequency its imag¬ 
inary part defines the elastic scattering rate of quasi¬ 
particles, T = -Im[E fl (0)] = n i v? mp Tr/(2V)'}2 ps A s (p), 
where A s (p) = — (l/n)lmg R (p,uj = 0) is the spectral 
function of each helicity band. In Fig.3a we plot the nu¬ 
merical self-consistent T as a function of the ratio n/no 
for different values of the SO coupling and we compare it 
with Boltzmann result r# = l/(2r). As expected Boltz¬ 
mann result is accurate for large n/no , where F <C Ep 
holds and one can approximate the spectral function as 
A s (p) — 8(Ep — Ep), and its accuracy increases with in¬ 
creasing Eq/Tq- On the contrary, as Ep approaches the 
band edge, the DOS singularity is smeared by disorder 
and finite-band effects cut-off the divergence of Boltz¬ 
mann result T b = To^Eq/E p following from (6). In 
the extreme diluted limit, (shaded region in Fig.3a) T as 
given by the SCBA vanishes and the diffusive approxi¬ 
mation breaks down [43]. 

Within linear response theory the conductivity is given 
by Kubo formula and, at T = 0, it is related to the on- 
shell current-current response function as follows [1] 

vdc = ^ {Pi R - Re[Hf*]) , (14) 

where P xx 4 is given by 

p if = y H p Tr [ix(p)G I ( P ,0)Jr(p)G M (p,0)] , 

(15) 

and the superscripts L,M = A,R indicate ad¬ 
vanced/retarded quantities. In the above equation j x ( p) 
and J X M { p) denote the bare and dressed currents and 
they are in general represented by 2 x 2 matrices in the 
helicity space. In particular, j x (p) = ev x . is proportional 
to the bare velocity (7), while J x has to be determined 
self-consistently [43] and it can be written as 



V±, M 5 aa > - ia LM s( 1 - 8 aa ,)t. 


(16) 


where V ps M = p/m + sa LM p denotes the dressed quasi¬ 
particle velocity. By comparing Eqs. (7) and (16) one 
sees that, as usual [1] in Kubo formalism the effects of 
scattering by impurities are encoded, via the vertex func¬ 
tion in the renormalization of the velocity. As we 

show below, under appropriate conditions, these effects 
are equivalently accounted for in Boltzmann language by 
the transport scattering times. 

The anomalous velocity a LM also plays an important 
role in the spin-Hall effect [46]. In this context it was 
shown that, although a RA = 0 in the HD regime [36, 47, 
49] and a RA ^ 0 in the DSO regime [48], in both regimes 
the spin-Hall conductivity vanishes. This result follows 
straightforwardly from the vanishing of d RA in the HD 
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regime while it can be proven by an explicit calculation 
in the DSO regime [48]. 

Using Eqs.(7) and (16), the current response function 
(15) can be cast as the sum of inter- and intra-band 
terms: P R X M = P R ™ a + P R ™, wh ere 

P-Lu, = (17) 

PLto r = E V , 9s (p, 0)9s' (p, 0) • (18) 

From a numerical self-consistent solution of the self¬ 
energy and vertex equations, we calculate the fully quan¬ 
tum dc conductivity Eq.(26). The results are shown in 
Fig. 3(b), where we plot the conductivity as a func¬ 
tion of the electronic density and of the SO coupling. 
Here we also plot Boltzmnann conductivity (blue lines) 
for two values of Eq, showing that Kubo results follow 
quite closely Boltzmann prediction. 

The equivalence between the two approaches can be 
proven in the the limit of vanishingly small broadening of 
the spectral functions, T <C E F , where we can discard[l] 
the RR term in Eq. (26). Indeed by neglecting also the 
inter-band contribution to P RA , relevant only at n — 
no where the spectral functions of the two chiral bands 
overlap, we can recast the conductivity as 

^ PS ^ ' ^ 5(Ef ~ ( 19 ) 

that, by direct comparison with (8), yields a*. ~ a R c 
provided that V RA ~ t^ s v ps /t. In the limit T <C Ep this 
relation is a straightforward consequence of the vanishing 
of a RA for Ep > E 0 and it can be easily proved for 
Ep < E 0 using a RA ~ a(l — Ep/E 0 ) [43]. This shows in 
particular that on the Fermi circles V ps /vp = r* r /r. 

The deviations between Boltzmann and Kubo results 
are better seen in Fig. 3c, where we compare the 
corresponding mobilities, respectively defined as fip = 
<T fc/( en ) an d Pt = cTd c /{en), and they can be ultimately 
ascribed to two factors. First, finite-band effects, that 
are mostly relevant for n < 0.3?ro and are responsible for 
the deviations of T from T#, shown in Fig. 3a, and for 
the relevance of RR terms [43] that in turn imply that 
for small densities fit tends to saturate, in contrast to 
fi b • Second, inter-band terms that are mostly relevant 
at n ~ no and give a smoothening of dependence of fit 
on n. This effect could be also captured by replacing 
the semiclassical Boltzmann equation with a fully quan¬ 
tum kinetic equation that includes also the off-diagonal 
components of the velocity operator (7) and of the non¬ 
equilibrium density-matrix in the helicity space. This 
allows one to account, in the presence of external fields, 
for the coherent superpositions of states with different 
helicities, as explained e.g. in Refs. [9, 10, 33, 50]. 

In conclusion, we have shown that in Rashba 2DEGs 
SO coupling entails an unconventional dc conductivity, 




Figure 3: (Color online) (a): Scattering rate in units of Fo 
as a function of n/no for different Eq. The dashed line 
shows Tb = l/2r (Eq. (6)). (b)Numerical dc conductivity 
(Eq. (26)) as a function of Eq/Tq and n for Fo = 0.5meV and 
m = 0.7 m e (as appropriate e.g. for LaA 103 /SrTiC >3 inter¬ 
faces [24]). The blue lines show Boltzmann conductivity for 
Eo = 50 and 100 Fo (c) Mobility fit (solid line) normalized 
to fit = e/{2hTom), compared with Boltzmann’s fiB (dashed 
line). Shaded areas in (a) and (c) denote the low-density re¬ 
gions beyond the diffusive approximation. 


strongly dependent on the density and on the SO- 
coupling strength. The experimental verification of our 
results requires the condition for diffusive transport 
(Ep T) to be fulfilled in the DSO regime (E 0 > Ep): 
this ultimately implies E 0 > Ep T. The conductivity 
anomalies in the DSO regime can then be accessible 
experimentally in relatively clean (T 0 ~ 1 meV) sam¬ 
ples of the strong-Rashba materials mentioned in the 
Introduction (see Refs. [3-32]), where Eq « 10 — 140 
meV and m ~ 0.2 — 0.7 m e , corresponding to electron 
densities no ~ 0.6 — 8 x 10 13 cm -2 [53]. Finally, we 
also remark that the large value of Eq in these systems 
guarantees that our zero-temperature results will provide 
a good description for real materials up to temperature 
scales kpT ~ E 0 . In addition, the occurrence of the 
DSO anomalies at relatively large densities n ~ no 
justifies also neglecting of electron-electron interactions, 
even though a full understanding of the conductivity 
anomalies in the diluted regime where interactions 
become relevant is certainly an interesting topic for 
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HELICITY EIGENSTATES BASIS 


To fix the notation, let us start by giving some details 
on the helicity eigenstates basis. The helicity operator 
is defined as S = [p x a\ z and its eigenstates, |ks), with 
s = ±1, satisfy the relation S'|k±) = ±|k±). A sim¬ 
ple calculation shows in particular that |k±) can be ex¬ 
pressed in terms of the standard spin eigenstates, |k f) 
and |k T)? as |k±) = (exp (—i$k)|k f) ± i|k \))/\/2 with 
0k = arctan(fcj / /fc ;c ). Consequently, the matrix 14 which 
implements the rotation from the spin to the helicity 
eigenstates basis has the form 


£4 = 


1 

71 




( 20 ) 


In the basis spanned by the states |ks) the total Hamil¬ 
tonian (Eq. (3) in main text) can be recast as follows: 


H — 7 Hk Ck + Emp(q)Ck +q ££k +q CkCk (21) 

k q.k 


where Ck = (ck+,Ck_) and = (c^ cj^) are spinor 
creation and annihilation operators, Hk is the Hamilto¬ 
nian of the clean Rashba model, i.e. H^ = diag(E ; (" — 
Eq,E// - E 0 ) with E± = ( k ±p 0 ) 2 /(2m), p 0 = ma and 
Emp(q) denotes the Fourier transform of the impurity po¬ 
tential, I4np(q) = 1/V e‘ q ' Rj, i)imp- Differently from 
the main text, in this Supplementary Material where not 
differently specified we use units h = e = 1. 


GREEN’S FUNCTION 


The Green’s function obeys the standard Dyson equa¬ 
tion: G _1 = (G°) 1 — £, where G° is the Green’s func¬ 
tion of the Rashba model in the absence of disorder and £ 
is the self-energy. Specifically, in the helicity eigenstates 
basis, G° has the form 


— (i^i — Ep + Ep) 1 S a /3 


where e; is a fermionic Matsubara frequency, and from 
now on we set the zero of the energy to —Eq. Thus, as in 
the main text, the dominant spin-orbit (DSO) regime is 
identified by Ep < E 0 . Within the self-consistent Born 
approximation (SCBA), E = E(p, ie n ) is determined by 
solving the following equation 


E(k, ie n ) = 


UiV 


imp 


UlU k G(p,ie n )UlU p ( 22 ) 



Figure 4: Wigwam diagrams which describe the self-energy 
within Born approximation. The solid line correspond to the 
dressed Green function, G, while the crosses indicate averag¬ 
ing over disorder [1, 2]. 


which corresponds the “wigwam diagram” depicted in 
Fig. 4 as described e.g. in Refs [1, 2]. Note that in the 
helicity basis to each impurity-scattering vertex, chang¬ 
ing the electron momentum from k to p, one has to as¬ 
sociate the spin rotation U^U^. Equation (22) admits 
a momentum- and spin- independent solution. Indeed, 
assuming [E(p, «e n )] Q/3 = E(ie n ), the momentum depen¬ 
dent part on the r.h.s. of this equation averages away. 
Within SCBA the Green’s function G(p, ie n ) is thus rep¬ 
resented by the following diagonal matrix in the helicity 
eigenstates basis: 

[G(P, fej)] aj 8 = (4; - E“ + E f - E(!ie,))" 1 ^. (23) 

By analytical continuation to real frequencies of Eq. (22) 
(see e.g. [2]) we obtain the following self-consistent equa¬ 
tions for the scattering rate T 

Tl ■ 7 2 F 

r = —Im[E*(0)] = E [\g+(p, 0)| 2 + \g«( P ,0)\ 2 } 

(24) 

and the retarded self-energy E fl (w) 

E 0?(p, - p) (25) 

where g±(p,u) = [cu — E± + Ep — Y, r (uj)] 1 is the 
Green’s function of each chiral eigenstate. To simulate 
a finite Brillouin zone, in Eq. (25) we introduced an 
upper momentum cut-off, p c . The latter is needed, in 
particular, to regularize the real part of the self-energy, 
Re[E^(w)], which would otherwise diverge logarithmi¬ 
cally at the band edge, see e.g. Ref. [3]. In these re¬ 
gards, we notice that, contrarily to what happens in 
standard half-filled systems where Re[E fl (u;)] is approx¬ 
imately ^-independent and it can be absorbed in a re¬ 
definition of the Fermi level, in the low-doping regime 
investigated here Re[£ fl (a;)] acquires a non-trivial fre¬ 
quency dependence. We thus need to calculate self- 
consistently both Re[E fl (o;)] and Im[£ fl (cu)]. Such self- 
consistent solution identifies the elastic scattering rate, 
r = — Im[E fl (0)] and the renormalized density of states 
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Figure 5: Structure of the real and imaginary part of the 
self-energy close to the lower band-edge for Eo = 40Fo. 



+ .. 


Figure 6: Ladder diagrams describing the conductivity within 
Born approximation. Solid lines and empty circles represent 
respectively G(p,iei ) and j x ( p). 


current response function in Matsubara frequencies: 
P xx (iei,i£i +n ) = ^ TV {G(p, isi)j x (p)G(p,i£i+„) ■ 

P 

■ J x (p,i£hi£i+n)j (27) 

where j x is the bare velocity operator introduced in the 
main text. We recall that in the helicity basis j x = ev x 
is represented by the following matrix: 


(DOS) N(E) = The electronic 

density at T = 0 is given by n = J E ^ N(E)dE. At 
low doping the impurities lead to a smearing of the van- 
Hove singularity in the DOS, which reflects in the be¬ 
havior of the scattering rate T, as described in the main 
text (Fig. 3a). In addition, the presence of impurities 
gives a shift, A e dg e , of the lower band edge, so that the 
lower “effective” band edge where n = 0 is identified by 
Ep = Ep + Aedge = 0-In Fig.5, as an example we show 
the structure of the real and imaginary parts of the self¬ 
energy as functions of the frequency, for frequencies close 
to the lower band-edge, located at Ep = —A 0 d ge - No¬ 
tice that these two quantities are connected by Kramers- 
Konig relations. 


DC CONDUCTIVITY FROM DIAGRAMMATIC 
PERTURBATION THEORY 

Current response function 

Within SCBA the static conductivity is given by 
Eq.(14) of the manuscript, that we report here for con¬ 
venience 


v x = — (To + acosf?p<T 2 + asin8 p a v (28) 
to 

where Co is the 2x2 identity matrix and cr, with i = x, y, z 
are the Pauli matrices. The renormalized charge current 
J x { p, i£i,i£ m ) satisfies the diagrammatic equation shown 
in Fig. 7 that in helicity space can be written as 

2 

J x (k,i£i,is m ) = j x { k) + 1 ™ P ^2 [ulU p G{p,i£i)- 

P 

J x (p,i£i,i£ m )G{p,i£ 

m )UlU k ] (29) 

Before coming to the solution of the above equation, 
we remark, that, as usual [2], the AR and RR response 
functions appearing in Eq.(26) correspond respectively 
to P xx (0 — iS, 0 + iS) and P xx (0 + iS, 0 + iS). 

As it can be easily verified, by symmetry arguments 
one finds that J x has the same matrix structure of the 
bare current (28), so that we can write: 

Px 

J x (k, iep ie m ) = —ap + a cos 8a z + a sin 9<j y (30) 

TO 

where a(i£pi £ m ) satisfies the following self-consistent 
equation: 


Vdc = ^ (Pi R - R e[P RR }) . (26) 

The derivation of the above equation in the absence of 
spin-orbit coupling is standard textbooks material (see 
e.g. Refs. [1, 2]) and, since it does not change in the 
presence of spin-orbit coupling, we do not review it here. 

By applying diagrammatic perturbation theory, one 
easily sees that the calculation of ad c implies the sum¬ 
mation of all ladder diagrams shown in Fig. 6. This 
in turn corresponds to calculate the following current- 


2 

T1 1). 

a = a + ' ^ mp ^ Ti[U^a y UpG(p, i£i) ■ 

P 

■Jx{p./i£l,i£m)G(p,i£ m )]■ (31) 

= > + D 15 

Figure 7: Diagrams describing renormalization of the charge 
current vertex. 
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Equation (31) can be solved explicitly to obtain to the 
following result for the renormalized anomalous vertex: 




a + cxq (is i, is m ) 
1 - A(iei,is m ) 


(32) 


where we introduced the quantities A(i£i,i £ m ) and 
a 0 (i£j,*e m ) given by: 


A(iei,i£ m ) = } j g s (k,i£i)g s '(k,i£ m ), (33) 

k ss' 

a 0 (i£i,i£ m ) = ' ,“ nP E -sg s {k,iei)g s {k,i£m). (34) 

4V ' m 

k s 


By replacing Eqs.(28) and (30) in Eq. (27), we arrive 
at the following expression for the correlation function 

P X X (i'-l , 771) ) 


F*XX (isl 5 ^ m ) 


P 0 (i£i,i£ m ) + 


m [(a + a)ao+aaA] 


(35) 

On the r.li.s. of the above equation the frequency depen¬ 
dences of a, «o and A, defined in Eqs.(31-34), is implied 
and we introduced the function Po(i£i,i£ m ), 


P{i {^£1 , i>£mri) 


1 x k 2 

■^z2—9s(k,i£i)g s (k,i£m), (36) 


2V *—' in 

k s 


that in the absence of Rashba coupling yields the only 
non-vanishing contribution to the conductivity. 

Once the analytical continuation is performed, Eq. (35) 
is equivalent to Eqs. (17-18) of the main text. This 
can be easily seen by inserting in Eq. (35) the explicit 
expression of A, «o and Pq, given above. One then finds 
that P R X M (0, 0) can be recast as 


pLM _ 

XX 



p 


p(a + a) 
m 



4 +««“') (d»f + 9 L -9 M ) 

m z J 

9 L _9 M ) + aa LM (g L + 9 M + 9 ^) 


From the above equation we see P xx is the sum of an 
inter- and intra-band contribution, i.e. 


pLM __ J_ 
xx 2V 


£{£[(^~)(^“) 


9s <?f 


+aa 


LM 


E 

s^s' 


9s 9 M 


’} 


= R 


LM 


P: 


LM 


(37) 


Eventually, since (^ + scc)(^ + sa LM ) = v ps ■ V ps M we 
arrive at Eq. (18) of the main text. 


Analytic approximations in the weak-disorder limit 

Starting from the above results in this section we de¬ 
rive approximate analytical expressions for the renormal¬ 
ized vertex and the conductivity. We assume that we are 


in the weak-disorder limit (WDL), where we can (i) ap¬ 
proximate the spectral functions with a delta, i.e. set 
A±{p, 0) = (T/tt)^ (p,0)| 2 = 6{E± - E f ); (ii) neglect 
the RR contributions. 

The WDL approximation for a RA can be derived start¬ 
ing from the analytic continuation of Eqs.(32-34). Per¬ 
forming the angular integral, A RA and ot RA can be then 
cast as 

ii ■ ii? r°° 

l P 2 [M 2 - \9 R (p, 0)| 2 ] dp, 

(38) 

A RA = J P 1S?(p, 0) + g R (p, 0)\ 2 dp. (39) 

The latter equation can be simplified using the self- 
consistent self-energy equation to obtain: 

11 11 ? fOO 

A RA ( 0,0) = 1/2 + jI pRe[g R (p, 0) g R (p, 0)] dp. 

(40) 

Approximating the spectral functions with a delta, as 
stated above, we immediately see that, except in a small 
density range around Ep ~ Eq, we can neglect the second 
term and the r.h.s. of Eq.(40), since there is no overlap 
between the two-chiral bands and we obtain A RA ( 0, 0) = 
1/2. In the same approximation, setting = E p — Ep, 
we can write a^ A (0,0) as follows 

ao A (°,°) = 8m ^ F j P 2 [<*(£p ) - S (Zp )] d Pi ( 41 ) 


which leads to 


a\ 


RA 


( 0 , 0 ) ~ - 


4 rrPvpT 


{p 2 _-p\) Ep>E 0 
(pl+pl) Ep < Eq 


with T 0 = niV? mp m/ 2. We recall that here, as in the 
main text, p± are the momenta on the inner and outer 
Fermi surface, so that their subscripts refer to the value 
of the transport helicity 77 , introduced in the main text 
in the context of Boltzmann transport. The sign change 
on the r.h.s. of Eq.(42) is thus due to the fact that for 
Ep > Eq the two contributions come from the two chiral 
bands in Eq. (41), while for Ep < E 0 only the E~ band 
contributes, with a two-folded Fermi surface. Using the 
explicit expression p±, i.e. 

pEp>E 0 _ mvF -pp 0 an( j p E F <E o _ p 0 || mVF . ( 43 ) 


with vf = ij2Ep/m, we eventually obtain: 


<*0^(0,0) 


—o; Ep > Eq 

— {Eq + E F )/pq Ep < Eq 


(44) 


Here we also used the WDL results for T, that can be 
derived from Eq. (24) in the WDL where \g R {p, 0)| 2 = 



10 


(7t/T) 6(E± - E f ) as 


yWDL _ f To _ E f > E 0 

\ Tq^/Eo/Ef = T 0 po/(mvF) E f < E 0 

(45) 

that coincides with the Boltzmann result from Eq. (6) of 
the main text. By replacing the result (44) into Eq. (32), 
along with A ra {Q, 0) = 1/2, we then obtain the estimate 
of a RA quoted in the main text, i.e. 


a flA (0,0) 


0 Ef > Eq 

a(l — Ef/Eq) Ef < Eq 


(46) 


Let us now discuss the analytic approximation of the 
conductivity. As discussed in main text in the WDL we 
can put 

pRA 

a dc ~ -jp (WDL) (47) 

Z7T 

where the intraband term P R // a coincides with the first 
term on the r.h.s. of Eq. (37). By using the result (46) 
for the anomalous vertex we can rewrite it as follows 

P&L = e f > E ° ( 48 ) 




RA = [p\ -P-- ma RA (p+ ~ p-j] 
intra 4?7i,r 


E f < Eq. 


(49) 

Using the expression of p±, T and a RA derived above, 
along with the expressions for the particle density in the 
WDL, i.e. 


_ f (m/7r)(E F + E 0 ) E f > E 0 

" l (p 2 o/n)^EFE~ 0 =n 0 ^EFE~ 0 E F < E 0 [ > 

Eq.s (48-49) lead to the final expression for the conduc¬ 
tivity quoted in Eq.s (l)-(2) of the main text. 

To conclude this section we would like to show that the 
inclusion of vertex corrections is crucial in both regimes. 
The “bare-bubble” conductivity a corresponding to 
the first diagram in Fig. 6, is given by the term P 0 de¬ 
fined in Eq. (36). It can be directly computed from Eq. 
(37) by replacing a with a, so that the renormalized ve¬ 
locity V ps is replaced by the bare one v ps . In the WDL 
we then easily obtain 

p u> A = FyE^.ls S ( E p- EF ) = ^(P++P~) ( 51 ) 

p,s 


Thus, using Eq.s (43), (45) and (50) into Ecp (47) one 
easily obtains that (restoring the charge e) 


Cbb 

Cbb 


e 2 (n — no/2) 


2 mV 


— ^Drude &n 


2 

e Tl 
4rnT 


—WDL 

a Drude 


E f > Po(52) 
E f > Eq (53) 


Figure 8: (Color online) (a) Different contributions to the 
conductivity as a function of the density for Eo = 40Fo , 
To = 0.5meU and m = 0.7m e . (b) Comparison of the total 
the AR mobilities. Parameters as in panel (a). 


where, as in the main text, dDrude = e 2 n/(2Tom), <r„ 0 = 
e 2 no/(4r 0 m) and = e 2 n/(2T WDL m). We there¬ 

fore see that the bare-bubble result is inadequate at all 
densities and chemical potentials. This also shows that 
even recovering the Drude conductivity at Ep > Eq is a 
non-trivial result, due to the crucial role of vertex cor¬ 
rections. Indeed, even in the regime Ep > Eq, where Ep 
slightly exceeds the Rashba energy, so that two Fermi sur¬ 
faces are clearly separated, any signature of the Rashba 
interaction disappears in the dc conductivity, which is 
given by the usual Drude formula. We also notice that 
vertex corrections tend to enhance the conductivity with 
respect to the bare-bubble result. As we shall discuss 
below, this is the result one usually expects within a 
Boltzmann picture, where backward and forward scat¬ 
tering processes contribute to the transport scattering 
time with different weights. 


Relevance of the RR contribution 

As we mentioned above, our analytical formulae are in 
principle valid only in the WDL, realized for T <C Ep ■ In 


2 
















this limit we did two approximations, we replaced the 
spectral functions of the chiral eigenstates with delta 
functions and neglected the RR part in Eq. (26). As 
the density decreases in the DSO regime these approxi¬ 
mations are not valid anymore and the quantum result 
for adc starts to deviates from the analytical WDL re¬ 
sult udso- In this regime the RR current response func¬ 
tion becomes as much relevant as the RA one, and it 
is important to include it in order to reproduce phys¬ 
ical results for the dc conductivity. In Figure 8(a) we 
plot separately the contributions of the RR and AR re¬ 
sponse function to the conductivity, defined respectively 
as a AR = P AR (0, 0)/(27t) and a RR = —P RR ( 0,0)/(27r). 
As one can see, for our parameters choice, at n/no < 0.5 
the RR contribution increases, while the AR one tends to 
saturates with decreasing density. Neglecting (jrr thus 
leads to the unphysical result of a finite conductivity for 
zero density. The vanishing of the conductivity as n —> 0 
is indeed guaranteed by the cancellation between a AR 
and <j rr . This fact is also evident in Fig. 8(b) where 
along with the total mobility gt we also plot the “AR” 
mobility, p AR defined as p AR = cr AR /n. We se that the 
inclusion of the RR contribution significantly modifies 
the structure of the mobility curves around the mini¬ 
mum. Eventually we note that the cancellation between 
RR and AR terms turns out to improve the agreement 
between adc and ajjso and it enlarges the range where 
the DSO formula can be used to describe transport. 


DC CONDUCTIVITY WITHIN SEMICLASSICAL 
BOLTZMANN APPROACH 


For a detailed discussion of Boltzmann equation we 
refer the reader to Ref. [4], here we only outline the most 
significant steps. Note that here we neglect the effect 
of the off-diagonal terms of the spin-density matrix in 
the helicity basis at a given p since as discussed in the 
main text they give only small corrections to dc charge 
transport (see e.g. Refs.[5, 6, 9, 10]). Let us call p PyS 
the distribution function for the quasiparticle eigenstates 
£p S , where s denotes in general the band index, which 
coincides in our case with the chiral index. The time 
derivative of p is determined by the collision integral: 

= ~ E Gp'Vbp.* - Pp',s'], (54) 

p',s' 


where Qp? s , is the scattering kernel from the state 
to the state £ p ' tS '- In the presence of an electric field E 
the l.h.s. of the above equations is given by: 


&Pp,s 

dt 


-eE • 


dpp.s 

dp 


„ - d PZ 

~ Vps de ps 


(55) 


in the last passage we replaced p with its equilibrium 
value p p q s = /(e ps ), since we are interested in the lin¬ 
ear response in E. In the relaxation-time approximation 
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we can express the time evolution of p via a transport 
scattering time r*£, so that: 


dpp.s 

dt 


Pp.s 


-pf q 


p,s 


(56) 


By combining Eqs. (55)-(56) we then have: 


Pp,s — P 


p.s 


+ eE • v 


Tl r 


dp eA 


PS 


P s 'ps 


de r 


(57) 


For a field in the x direction the current can then be 
written as: 



. (58) 

where we used the fact that there is no current in the 
equilibrium state. At T = 0 we can put <9p p q s / de ps = 
df(£ P s)/de ps = 6 (e ps -/x). Thus, by identifying e ps -p, = 
Ep — Ep we arrive at Eq.(8) of the main text 

^ = Y E \vps\ 2 rl r s 8{E ps - E F ) (59) 

ps 

A set of equations for the transport scattering times 
can be derived by substitution of Eq. (55) into Eq. (54), 
once that one uses the Ansatz (56). By doing so, since 
the scattering kernel Q conserves the energy, and the 
equilibrium function p p q s does not depend on the chiral 
index but only on the energy, one is left with: 


E O ps \r tr v ’ 

^pV L'ps l P s Vs u p 


ps 

Je; 


w ps 


-E 


’p's' 1 


(60) 


p's' 


where we introduced the quasiparticle scattering time 


r 





ps 

p's'* 


p's' 


(61) 


By using the fact that only the component of v p > s > in the 
direction of v ps survives after momentum integration, Eq. 
(60) finally reduces to Eq.(9) of the main text 





p's'" 


^p's' 


^ps 


p's' 


(62) 


Notice that Eq. (62) differs from the one proposed e.g. 
in Ref. [7], where the band-dependence of the transport 
scattering times on the r.h.s. of Eq. (62) has been over¬ 
looked, leading to decoupled equations for the T* r s . Here 
instead the set of coupled equations (62) is analogous to 
the self-consistence equations (29) introduced above for 
the renormalized current in the quantum language. This 
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analogy can be exploited further by the identification of 
the renormalized Boltzmann current as 


J 


B 

ps 


ps 


ps 


( Es p 


(63) 


As already discussed in Ref. [8] for the case Ep > E 0 , 
both the quantum and the Boltzmann approaches lead to 
the same renormalized currents. This result also holds 
in the DSO regime Ep < E 0 , as one can see from the 
explicit solution for the r’s derived below. 


Collision integral 


can rewrite the above equation as: 


ps 

r 


= i 


L imp \ ^ 


V 


ip ■ V') VpsVp's'S(E F - E 8 p ,)t p 


( 68 ) 

where we set r (Ep) = r and ij ps = s(v ps ■ p) = ±1. 
Performing the angular integral p' and changing variables 
from p,s to E,rj we eventually recover equation (10) of 
the main text: 



1 

AtqIJIVp 


v' 


(69) 


By using Fermi Golden Rule, the scattering rate from 
the state |p,s) to the state |p',s) can be written as: 

r ) , jr 

QVs' = y\(Ps\V imp \p's')\ 2 8(e ps - £pV ) (64) 

Using the explicit expression of the helicity eigenstates 
in plane waves, we can rewrite the above equation as 
follows: 


QPf = — 

^P s' p2 


y'dre i (P-P' )r U imp (r)lUPP' 


6(e;-e;,) 


(65) 

Here the matrix W pp = Up,U p comes from the scalar 
product of the helicity eigenvectors. Within our ap¬ 
proximations, (self-averaging delta-correlated disorder 
and Born scattering), we can write (V^mp(r)V^mp(r , )) — 
riivf mp 5(r — r') and we can recast the above equation as 


Q 


ps 

p's' 


2tt 2 

~y~ n i v imp 


W P 


pp 


he; 



( 66 ) 


where 


lr p P '|2 = 1 + sign(sa / )cos(gp 



so that one recovers Eq. (5) of the main text. 


Solution of Boltzmann equations 


Using the explicit expression of Q p s s , Eq.(62) reads 


(e; 


j imp 

V 


^(1 + ss'p-p') 

p's' 


■5{e;-eM s ,. 


(67) 


Now recalling that v ps = v ps p and that |u ps | = v F = 
y/2rnEp for E; = Ep , for states at the Fermi level we 


where tq denotes as usual the quasiparticle scattering 
time in the absence of spin-orbit, To = l/(mnjif mp ), p v = 
p v (Ep) indicate the two Fermi momenta introduced in 
Eqs.(43) and we set r* r = T p (E F ). 

To solve this equation it is useful to note that t/tq = 
mvp/pp where p F = 1/2 J2pPv that allows us to recast 
Eq. 69 as: 

= 1 + \ Y, VV'PpTp■ (70) 

n' 

with r* r = r‘ r /r and p v = p v /p F . As one can easily 
check, the solution of this equation reads: f* r = p v that 
coincides with Eq. (12) of the main text. 
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